Loss in hybrid qubit-bus couplings and gates 
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We provide a characterization and analysis of the effects of dissipation on oscillator assisted 
(qubus) quantum gates. The effects can be understood and minimized by looking at the dynamics of 
the signal coherence and its entanglement with the continuous variable probe. Adding loss in between 
successive interactions we obtain the effective quantum operations, providing a novel approach to 
loss analysis in such hybrid settings. We find that in the presence of moderate dissipation the gate 
can operate with a high fidelity. We also show how a simple iteration scheme leads to independent 
single qubit dephasing, while retaining the conditional phase operation regardless of the amount of 
loss incurred by the probe. 
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I. INTRODUCTION 

Quantum information processing can potentially pro- 
vide a considerable speedup over classical information 
processing for certain problems |l[ along with the ability 
to efficiently simulate physical systems that cannot be 
done classically @, S|- I n view of this, much work has 
been done on finding a viable physical implementation 
of a quantum computer and many different realizations 
have been proposed [4|, |5|, [fj, [D, |8|, some of which ex- 
plored experimentally on a small scale [fj, LLCj, [ll|, Il2j |. 
Most large scale architectures rely on protocols enabling 
the transport of quantum systems or distributed schemes 
to perform logical gates on isolated qubits. Such schemes 
make use of 'flying' qubits 



quantum bus concepts 




14J or more generally of 
Decoherence effects on 
the bus are crucial and being able to overcome them will 
bring us one step closer to true scalability. 

Most of the results in quantum information theory 
were developed in a discrete setting, making use of qubits. 
However even though single qubit operations are not 
so much of an issue, the level of control required to 
physically implement entangling gates between individ- 
ual qubits and measure them within realistic coherence 
times is tremendous, limiting the experimental realiza- 
tions. The initial theoretical proposals were rapidly 
adapted to a continuous variable (CV) setting, where CV 
quantum information processing was shown to be possi- 
ble [17|, [l8|] . CV implementations may be more accessible 
in some respects, with simple measurements and entan- 
gling operations. Despite these advantages, this frame- 
work is limited by the nonlinearity available experimen- 
tally, making single system operations difficult. Combin- 
ing the exactness of discrete variables and the robustness 
of continuous variables is therefore a judicious route to 
take. The term 'hybrid' was first coined by Lloyd [l9( to 
describe quantum information processes featuring both 
discrete and continuous variables. The ability to switch 



on and off particular interaction Hamiltonians enables 
one to simulate interactions and quantum logical gates 
on discrete systems. Multi-qubit extensions were rapidly 
undertaken in generalizations 2fJ, |21 1 . The initial obser- 
vation came at the same time as Milburn's proposal to 
simulate interactions between trapped ions by coupling 
them to a common vibrational mode [22J . constituting 
a direct physical realization of a hybrid quantum com- 
puter. This approach was used to entangle up to four 
ions experimentally 10]. In recent years other physical 
implementations have been explored making use of mat- 
ter qubits [la] . A hybrid quantum computer could poten- 
tially be more versatile than its strictly discrete counter- 
part, providing simple algorithms to compute eigenvec- 
tors and eigenvalues [l9( or implementing Grover's search 
algorithm in a direct way [20j. 

A CV bus (qubus) can also be used to achieve quantum 
non-demolition (Q ND) measurements and parity gates 
23, [24|, 12a, l2a, [23] , which for example can be combined 
with an ancilla qubit to implement a near determinis- 
tic CNOT gate [28j , In these cases the limited strength 
of the nonlinearity is compensated for by an increase in 
the bus amplitude, leading to entangling gates based on 
homodyne measurements functioning at greater success 
probabilities than in single photon mediated applications. 
This is particularly relevant to the generation of cluster 
states [27J, I29J an d entanglement distribution for repeater 
applications [30. l3l|. 

The physical circumstances in which one can envisage 
a hybrid coupling between a CV and discrete variable 
system have been extensively investigated. The Jaynes- 
Cummings model [321 ] is very successful at accounting 
for the interaction of radiation with an atom in a cavity 
quantum electro-dynamics (CQED) setting. Based on 
this, one possible realization for the qubus scheme is that 
of an atomic qubit interacting dispersively with a cavity 
mode. However superconducting charge qubits are also 
a good candidate for a physical realization la, |33j , for 



which a dispersive coupling with a microwave bus mode 
has already been experimentally demonstrated 3J, I35J . 
Decoherence effects during such an interaction have been 
explored in the past, for example in the case of a two-level 
atom interacting dispersively with an optical mode in a 



dissipative cavity [31| 



361 1371 ] . This dispersive interaction 
forms the basis for many qubus schemes and dissipation 
effects during the interaction on an entangling gate be- 
tween two qubits can potentially be overcome [38J ]. A 
symmctrization technique to develop resilience to both 
dissipation and thermal fluctuations was also proposed 
[391 ] . There the authors took advantage of the invariance 
under time reversal of the action of the gate, noticing that 
the combination of an interaction sequence with its time 
reversed version canceled out dissipation effects. Even 
though they are effective, these methods focus on a pair 
of qubits within the same cavity or trap and decoher- 
ence due to inter-cavity communication remains to be 
addressed. In addition to this they propose interaction 
sequences which have to be iterated many times before 
any significant improvement in gate fidelity can be ap- 
preciated. 

The cross-Kerr effect is also well known and used in 
several different areas of optics. In this context, deco- 
herence effects have been approximated along with their 
impact on optical quantum information processing using 
40 . Ul|. It was shown that dissipa- 



weak nonlinearities 
tion effects in a two-qubit parity gate can be minimized 
if one can implement a photon number measurement on 
the bus. In the present work we extend previous results 
by solving exactly the master equation during the inter- 
action for arbitrary input states and observing a trade 
off between the entanglement generated and the required 
precision in interaction time. We also provide the result- 
ing quantum operations incurred by qubits involved in a 
CZ gate assuming sequential interactions. These oper- 
ations are critical when quantum error correction tech- 
niques have to be introduced as the system is scaled up. 

We note here that direct decoherence effects incured 
by the qubits, independently from the bus, are being ne- 
glected. In so doing we obtain a simplified picture, which 
can become a good approximation if bus dissipation is 
the primary source of noise. The eventual optimization 
of a full architecture will nonetheless need to factor in 
additional sources of noise. 

The paper is structured as follows. In the first section 
of this paper we derive closed expressions for the effects 
of dissipation during cross-phase and dispersive Jaynes- 
Cummings interactions. We follow the entanglement and 
coherence dynamics of a qubit and the continuous vari- 
able as they interact with each other. In section II we 
carry our attention over to the hybrid gates themselves 
starting with the conditional displacement gate followed 
by the CZ gate in section III. By adding dissipation be- 
tween interactions we obtain the quantum operations un- 
dergone by the qubits and discuss gate fidelity. Finally 



we provide a simple iteration scheme to simplify the oper- 
ation down to a perfect CZ gate followed by independent 
single qubit dephasing. 



II. INTERACTION LOSS ANALYSIS 

As explained in the introduction we focus here on two 
types of interactions. The first consists in the cross phase 
modulation undergone by probe and signal modes as they 
travel through a Kerr medium. The second type of in- 
teraction we consider takes place between an off-resonant 
optical probe mode and a single atom in a cavity. Their 
effects are described by interaction Hamiltonians of the 
form 
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where a(a^) are the annihilation (creation) operators act- 
ing on the probe mode. In the case of a cross-Kerr in- 
teraction, A = b% the number operator acting on the 
signal mode (A|n) = A n |n) = n\n), the states \n) cor- 
responding to the Foch state basis) and \ IS propor- 
tional to the third order nonlinear susceptibility of the 
medium. For the CQED setting, the interaction is de- 
scribed by the dispersive limit of the Jaynes-Cummings 
Hamiltonian 42j in which A = Z where Z is the Pauli 
operator acting on the atomic qubit (Z = |0)(0| — |1)(1|, 
A|n) = A„|n) = ( — l)"|n)) and x is the atom-light cou- 
pling strength [27j. In both the cross phase and the 
CQED settings we have a hybrid interaction between the 
continuous quadrature variables of the probe field and 
the discrete degrees of freedom of the subsystem. Initi- 
ating the probe in a coherent state \a) and applying the 
interaction for a time t yields 
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FIG. 1: Loss in the probe mode during the coupling between 
the discreet system \ij;) and the continuous variable prepared 
in the coherent state |a). 

Given these interaction Hamiltonians, we can evaluate 
the effects of dissipation in the probe mode during the 



interaction (see Fig . 1) by solving the optical Linblad 
master equation [43j 

-&■ = -\ Wmt, p(t)} +-f(2ap(t)a t - a)ap{t) - p{t)a)a) 

(3) 
where we have assumed a zero temperature bath (a good 
approximation in the visible light regime) . The damping 
factor 7 quantifies the dissipation in the probe mode. 
If we consider a general input density matrix element 
|n)(m| ® | a) (a | in which the probe and the signal are 
disentangled, we can find the equation of motion for this 
particular element by looking at the operator p nm {t) — 
(n\p(t)\m). Due to the disentangled form of the initial 
state we have p nm (0) = \a)(a\ for all n and m. The 
equation of motion for each element is given by 

— ^t — = ix^na) ap nm (t) + -ix\ m p nm (t)a^a 
at 

+"f{2ap nm {t)a) - a) ap nm (t) - p nm (t)a)a). (4) 

Following the method used in 36|, we use the super- 
operators M(-) = a^a(-), V{-) = {-)a)a and J{-) = a(-)(tf 
to rewrite the above equation as 



dt 



{ixiKM - \ m V) + 7(2.7 -M-V)} Pnm {t) 
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(5) 



The formal solution to (5) is then p nm (t) = e Cnmt p nm (0). 
The super-operators realize an algebra obeying the com- 
mutation relations [J, M\ = [J, V\ = J for which de- 
composition theorems have been derived [44l |. leading to 



exp[£„ m i] = exp 
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x exp [(iXnX - i)Mt\ 
x exp [(— iX m x - l)Vt\ . 
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Now applying this result to our initial clement p nm (0) = 
|a)(a| we obtain: 



Pnm(t) = exp[-|a| 2 {l - e 



1 _ P (-27+i(A n -A m )x)t 



x \ae 



(-7-MA„x)t\ / ae (-~(+ iX ™x)t I _ 



A m )x/2 7 



(7) 



The coefficient derived above is the closed expression for 
the 'coherence parameter' given in [40J. If we denote 
this coefficient by £„ m then we have £„ m = Cmn an d so 
I Cum I = I Cmn I as well as Cnn = 1- We also quickly notice 
from (7) that this coherence parameter does not tend to 
as t tends to infinity for a fixed a. We have 



|S,nm|i — >oo CXp 



(A„ — A m ) z 
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(8) 



One way of understanding this is that the probe under- 
goes loss as it couples to the signal, thus reducing the co- 
herence in the signal. But eventually the probe returns to 
the vacuum state (the time it takes depends on the initial 
amplitude and the damping factor), disentangling itself 
and leaving some coherence in the signal. However the 
larger the amplitude, the larger the effective interaction 
time, the less coherence remains in the signal. 

As time progresses, the process can be viewed as a 
unitary operation between the signal and the damped 
probe in addition to a dephasing effect on the signal. To 
view this more clearly, let us write the output density 
matrix using A = Z, 9 = \t and defining z n = (—1)", 



P(t) = ^2 C nniCnm\n)(m\ 

n,m— 0,1 

<8> \ae-^ t+lz " e }(ae-^ t+lz '" \ 
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Writing ( n 
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the factor e Rc ^™ m l character- 



izes the decoherence and takes the form e - £ ( 1 - 2 n z m) 
(see the Appendix for the full expression). Applying 
this type of operation to a qubit density matrix p — 
Cnm | n)(m | yields directly the phase flip channel 



,,m=0,l L ' nm 



e Re[/„ m ] p 
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The additional phase e Im ^ nm ] acquired in the process (see 
Appendix) is known and can be corrected for if needed; 
it is not an intrinsic source of noise. 

One expects this issue of coherence to be intimately 
linked to the entanglement shared between the signal and 
the probe systems. In order to observe the dynamics 
of this entanglement we restrict the signal to being a 
qubit and will continue to use A = Z. We also take a 
real for simplicity. Equation (7) provides us with a time 
dependent density matrix and having our input signal in 
the state (|0) + |1))/V2, it reads 



P(t) 



1 



{|0)(0|®|a }(ao| + Coi|0)(l|<g>|ao)(ai 



Cio|l>(0| ® |ai)(a | + |1)(1| ® |ai)(ai|}, (11) 



with ao = ae < - _7+lx - 1 ' and a.\ = ae'" 7- ^''. The entan- 
glement being invariant under local unitary operations 
we allow ourselves for simplicity to apply the conditional 
phase |0)(0|+e lIm [ aoQ! i]|l)(l| on the qubit. Then we rede- 
fine the bus probe states as |ao) and \a[) = e llra -[ a o a i\ | ai } 
so that the overlap between the two is real: (aolai) = 
|(cto|ai)|. This allows us to express them in an orthogo- 
nal basis {\x), \y)} as [46| 



\a ) = a\x) +b\y), 

K) = <A x ) -~ b \y)- 



(12) 



Taking a and b real without loss of generality, normaliza- 
tion leads to 



1 



b = 



1 



(13) 



(i cos(2 X t)) ^ this point we can write 
our locally equivalent density matrix in the orthonormal 
basis {10)11), \0)\y), \l)\x), \l)\y)} as follows 
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We have now managed to express the qubit and con- 
tinuous variable composite state in the form of a two 
qubit state. Given the resulting two-qubit density ma- 
trix, there are several entanglement measures to choose 
from, including the logarithmic negativity and the rela- 
tive entropy of entanglement 47J . Here we will work with 
the concurrence as defined by Wooters 48|, |49( , which we 
plot as a function of the scaled time %t for particular 
choices of parameters a and "f/x- As the qubit and the 
field initially start in a product state, and eventually for 
large times should return to a product state when the 
probe field doesn't contain anymore photons, we would 
expect the entanglement to peak at some point in time. 
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FIG. 2: Plots of the concurrence C(p) (solid) and the von Neu- 
mann entropy S(p) = — tr(plogp) (dashed) of the combined 
state of the continuous variable mode and qubit as a function 
of the scaled time xt- ( a ) From left to right a=200,100 and 
50 with a fixed ratio of damping rate to nonlinearity 7/x — 1- 
(b) The amplitude a is fixed to 100 and from lowest to highest 
peaking curves 7/x=l> 7 an d 21. 

This is verified in the plots of Fig. 2. In Fig. 2(a) we 
can observe how the amplitude of the field a affects the 
behavior of entanglement in time. As expected the larger 
a is, the larger the maximum entanglement. This is sim- 
ply explained by the fact that the separation between the 
possible states of the field in phase space increases with 
a, for a same interaction strength, thus making them 
more distinguishable. For large a the maximum concur- 
rence tends naturally to 1, however the peaking of the 



entanglement also becomes sharper. It is quickly gen- 
erated, but also quickly destroyed, as illustrated by the 
von Neumann entropy S(p) characterizing the decoher- 
ence which is a function of the squared distance in phase 
space between the two field states. Fig. 2(b) shows us 
how the maximum achievable entanglement depends on 
the ratio 7/Xj as did the limit of the coherence parameter 
in the previous section. The larger the relative damping 
7, the lower and the quicker the entanglement peaks in 
time. In both plots S(p) tends to 1, meaning the qubit 
is left in a maximally mixed state, disentangled from the 
probe. However this is not always the case and in general 
the smaller the ratio x/l i s > the lower the final entropy 
of the qubit becomes. 

In view of a QND measurement on a single qubit 
5fJ, [51 1, only the entanglement with the probe needs to 
be taken into account, as decoherence in the process will 
not affect the measurement statistics. However, when 
the application becomes cat state generation or multi- 
qubit gates 40j, |4lJ , decoherence becomes a crucial issue 
which has to be weighted against the entanglement. In 
such applications one wishes to produce coherent super- 
positions of single or multiple quantum systems. Thus 
it is important to view the behavior of the entropy of 
the combined state at the time at which the entangle- 
ment is maximized. This behavior is illustrated in Fig. 
3(a), showing the expected limiting behavior of the max- 
imum entanglement as a increases. In Fig. 3(b) the 
corresponding entropy of the combined outgoing state is 
seen to decrease asymptotically for all choices of the pa- 
rameter 7/x- In consequence one can simply reduce the 
amount of decoherence by increasing the strength of the 
probe. This is in part due to the fact that the interaction 
time becomes shorter, reducing the effective decoherence 
time. 

The success of such an approach to minimize the deco- 
herence will then depend on the loss incurred in between 
interactions. The reason being that the larger the ampli- 
tude a is, the larger the amount of dephasing incured by 
the qubits coupled to the probe mode during these time 
intervals will be. This will become clear in the next sec- 
tion. In consequence we observe a similar trade-off of as 
that encountered in schemes such as the hybrid quantum 
repeater proposed in [311 ] - If the transit time and con- 
ditions are appropriate, the approach is indeed effective. 
For example taking 7/x = 5 and a reasonable amplitude 
a = 10 4 we obtain a maximum concurrence of 0.998 for 
a von Neumann entropy of 10 -2 . 

So the higher the entanglement we want to measure 
or couple out of the cavity if we are dealing with cavity 
QED systems, the larger the probe amplitude and the 
more precise the timing of the interaction will have to 
be. Clearly, these issues of coherence and entanglement 
will have to be combined in order to optimize quantum 
gates in which different qubits interact with the same 
probe mode. 
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FIG. 3: (a) The maximum concurrence as a function of the 
amplitude a of the probe with y/x = 1, 3, 5, 10 and 15 
from top to bottom, (b) The von Neumann entropy of the 
combined state at the entanglement peaking time with the 
same values for y/x m decreasing order from top to bottom. 



III. THE CONDITIONAL DISPLACEMENT 



271 . l30l . l38j rely directly 



Many qubus gate proposals [If 
on the dispersive Jaynes-Cummings interaction. In those 
schemes different qubits in separate cavities interact se- 
quentially with the same probe mode and subsequent ho- 
modyne measurements can project the qubits to entan- 
gled states in a heralded fashion. Alternatively photon 
number resolved detection can be used, potentially lead- 
ing to deterministic gates. Dissipation in the probe mode 
during such gates will reduce the fidelity of the post- 
selected entangled states and has been investigated in 
some detail [31| . An interesting trade-off arises between 
the distinguishability of the measurement outcomes and 
the amount of dephasing incured by the qubits. 

More in line with the initial hybrid proposals [2fJ, 



21 



22j, measurement- free quantum gates can be imple- 
mented with the use of what we will refer to as the con- 
ditional displacement interaction. A qubit conditionally 
displaces a continuous variable bus with the following 
operation: D(j3Z) = exp[(/w — j3*a)Z] where j3 is the 
amount by which the field is displaced. Thus we can 
see that after the operation, the position of the probe in 
phase space is dependent on the state of the qubit (they 
have become entangled). 

There are different ways of arriving at such a coupling 
between a qubit and a continuous variable. The physical 
system may for example be naturally described by an in- 
teraction Hamiltonian directly producing the conditional 
displacement as defined above. A flux qubit in a cavity 
for example could displace the cavity mode in such a way 
|52j. In cavity QED systems, a particularly shaped pulse 
entering and leaking out of the cavity could also provide 
us with a very similar interaction Hamiltonian 53] . How- 
ever, based on the dispersive Jaynes-Cummings interac- 
tion, it was shown lately that a conditional displacement 
as defined above can be exactly simulated by using un- 



conditional displacements 5_4j. Defining the conditional 
rotation operation encountered in the previous section as 
R(9Z) = cxp[i8a^aZ], a possible sequence of interactions 
is the following 



D(acosd)R(-6Z)D(-2a)R(9Z)D(acos6) = D(2iasin6Z), 

(15) 
with a real (see Fig. 4) . Surprisingly all the phases cancel 
out to give us an exact conditional displacement, inde- 
pendent of the initial probe position in phase space. Mak- 
ing a pure imaginary will lead to a displacement in the 
orthogonal (real) direction. Consecutive displacements 
induce phases as D{(3)D{a) — e lIm ( Q ^D{a+(3) and con- 
sequently L>t(a)TJ»t( / 3) = e - ilm ( a *VD^(a+{3) as D^(S) = 
D(—6). Displacements and rotations acting on the den- 
sity matrix we will denote as T>{a)p = D(a)pD' f (a) and 

n(o) P = R(o)pitf(o). 
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FIG. 4: (a) The interaction sequence for the the simulation 
of the conditional displacement gate, (b) The two possible 
trajectories of the probe in phase space, the upper and lower 
paths corresponding to the qubit states |0) and |1) respec- 
tively. Here as in the text, we assume a to be real. The 'R' 
and T exponents denote the real and imaginary parts of the 
probe amplitude S at each time step. 

In order to characterize loss in an interaction sequence, 
we will introduce dissipation in the probe mode between 
and during each interaction. Dissipation during the inter- 
actions is due to loss in the nonlinear material or cavity 
system used to to mediate the interaction, whereas dis- 
sipation in between each interaction is due to fiber loss, 
mode mismatch and other effects. Consequently a differ- 
ent loss parameter should arise. However, both types of 
loss result in a dephasing of the qubit (s) and their effects 
can thus be combined into a single loss parameter. We 



will also assume the amount of loss is the same in each 
segment. Dissipation affects a coherent state matrix ele- 
ment as £|a)(/3| -> (P\a) r >\d){$\ with t) = 1- e" 2 ^ $j%. 
We fix the amount of loss between each interaction I = yt 
and represent the attenuated coherent state by d = ae 
such that the number of dots will determine the number 
of attenuations i.e. a — ae~ 21 . The quantum opera- 
tions on the qubits will be obtained by calculating these 
state-dependent overlaps. 

Now let us consider the effects of dissipation in the 
whole interaction sequence (15). For generality we will 
keep the amplitudes of the three displacements as free 
real variables aj., a 2 and CK3. The first step in the se- 
quence is a displacement so taking our probe initially in 
the vacuum state we have 



V{a x )p = J2 cw|a)(a'|<8>I>(ai)|0}<0|I>t(ai) 

a, a/— 0,1 

= ^2 c aa '\a){a'\ <8|ai)(ai|, (16) 



with T = sin#(dia2 + d'id2 + (a 2 + a 2 + d|)cos#) and 
S = sin 0(df + a\ + a?,). Here as in the evaluation of the 
coherence parameter (7) in section I, the exponent can be 
separated into real and imaginary parts. The former is 
characterized by S, representing the amount of dephas- 
ing incurred by the qubit, which can be decomposed into 
a phase flip channel as in (10). The latter constitutes the 
known phase acquired by the qubit in the process, charac- 
terized by T. Interestingly, this overall conditional phase 
can be tuned at will by adapting the amplitudes of the 
displacements, leading to the exact simulation of a con- 
ditional displacement. In other words, in this sequence 
we can limit the effects of dissipation to a dephasing on 
the qubit, producing nonetheless the correct combined 
output state. Previously also, loss during the interac- 
tion had led to single qubit dephasing effects which can 
be factored in here. We now move on to examine a two 
qubit gate in the presence of probe loss. 



where {|a),a = 0,1} represents the basis states of the 
qubit and the state of the probe mode is kept to the 
right. Loss in the probe mode at this point will decrease 
the amplitude of the coherent state to |di), without af- 
fecting the qubit. Proceeding to the second step which 
is a conditional rotation we have 



K{6Z).C.V{a 1 )p= Y, coa'|o)(o'|®|die w *»)(di 



e i6z a , 



(17) 



Now we introduce loss in the probe mode, leading to 
£.K(eZ).V(a 1 ) P = Yl W|a)(a'|0|d ie ^)<d ie ^'| 

a, a'— 0,1 

x \a\e a \a\e ) '. (18) 

Continuing in this fashion and completing the sequence 
C a p = V{a 3 ).£.K{-9Z).£.V{a 2 ).£.K{9Z).£.V(a 1 )p 

= Y c aa'\a)(a'\ <8> I a\ + a 2 e~ l0Za + a 3 )(z a -»■ z a >\ 

a, a'— 0,1 

xexp[isin6>(d 2 a3 - a 1 a 2 ){z a - z a >)] X (£1663)'', (19) 

with the £'s representing the three loss terms (overlaps). 
The notation z a — > z a > means the contents of the bra 
are the same as in the previous ket replacing z a with 
z a i . In order to simulate a conditional displacement (with 
dephasing on the qubit), we require the state of the probe 
mode to be of the form \jz a )(jz a i\. This is achieved by 
setting a'i + o^cos^ + a% = 0. Combining the loss terms 
we obtain 



666 = exp[-5(l - z a Za')] x exp[zT(z a - z a >)], (20) 



IV. THE CZ GATE 



Based on four of these conditional displacements in- 
duced by a pair of qubits on the same probe mode, 
a gate locally equivalent to the CZ gate can be built. 
An interaction sequence leading to this unitary oper- 
ation is U = D{-f3 b Z b )D{-[3 a Z a )D{f3 b Zb)D{p a Z a ) = 

e 2iIm{P*0 h }Z a Z h with fit-fa = i7T /g where Za and Zb act 

on qubits a and b respectively. Here we use the same 
method as before, introducing dissipation between each 
interaction, as illustrated in Fig. 5. We will assume the 
probe starts off in the vacuum state, but the same result 
is obtained if we initialize it in any coherent state. The 
initial density matrix reads 



Pab = 



= J2 c aa , w \ab){a'b'\ ®\0){0\. 



(21) 



a,a',6,&'— 0,1 



At this point it is worth noting that even though the 
probe undergoes amplitude damping, it can nonetheless 
be perfectly disentangled from the qubits coupled to it 
whatever the amplitude of the coherent state it starts off 
in (it may not be the case for other optical states). This 
is done by tuning the second (opposite) conditional dis- 
placement, in function of the known loss parameter I into 
which loss during the interaction has also been factored. 
Thus the amplitude of the second conditional displace- 
ment will be reduced by a factor of e~ 2 ', the damping 
undergone by the probe since the last coupling. Resolv- 
ing the whole gate sequence (see Fig. 5) and choosing (3 a 
and p b to be real we obtain 




FIG. 5: The interaction sequence for the CZ gate with 
loss. Between each one of the four conditional displacements 
(1,3,5,7) the probe undergoes dissipation (2,4,6). 
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xexp[i(/3 Q /3 b + (3a/3 b )(zaZ b - z a ' *%')]• 



(22) 



First we notice the geometrical phase represented by the 
last term. This is precisely the form of a two-qubit con- 
ditional phase having been applied to the density matrix. 
Ignoring the other terms, if we can set $ a Pb + (3 a $ b = 7r/4 
then we have simulated a CZ gate. However the dephas- 
ing effects are included in the three £ overlaps. The first 
and third lead to single qubit dephasing 



6 = (PaZa'IPaZa) = exp[-/3„(l - Z a Z a >)}, 

6 = (i$ b z b '\i/3 b z b ) = exp[-/?j(l - zbZb')]> (23) 

while the second overlap corresponds to loss in the probe 
mode when it holds information on both qubits 

6 = (PaZ a > + i(i b z v \(i a z a + i(3 b z b ) 
= exp[-0l + 01) + (/3 a z a + if3 b z b )((3 a z a , - i[3 b z b ,)}. 

(24) 

In order to be able to express the resulting quantum op- 
eration in a closed form, we first arrange the terms in 
the exponential such that when writing the expansion 
we obtain an accessible closed algebra. Symmetrizing 
the terms we have 



£,' = exp[—Xo + XlZ a Z a '+X2ZbZb>+X3(z a + iZb)(z a '—iZv)], 

. . ( 25 ) 

With X = 7]0a + /?&)> X l = VPaifla ~ Pb) , %2 = T}0b(Pb ~ 

(3 a ) and £3 = r]f3 a f3 b . Focusing first on the x$ term, we 
can write out the expansion of that particular exponential 



FIG. 6: The trajectories of the probe state in phase space. 
All four paths start and end at the same amplitudes a and 
a respectively, insuring that the probe disentangles from the 
qubits. 



acting on the density matrix as 

2- n\ ■ (26) 

n=0 

Now we define the two-qubit operators 

J = Z a + iZ b , 

K = Z a Z b . (27) 

Then grouping the terms in the expansion we obtain the 
following unnormalized operation 



a x 3 {z a +iz b )(z a i ~iz b t) 



,] p = c+p+C-KpK+s-JpJ^+s+J^pJ, 

(28) 



with 



c± = (cosh2:E3 ± cos2a;3)/2, 
S± = (sinh2a;3 ± sin2x3)/4. 



(29) 



Now we have to factor in the other terms xq, x\ and xi 

>2 ' 



of £2 ■ Further identifying 



cq = coshaqcosha^, e% = coshaqsinha^, 

e2 = sinhxiCOshx2, e 3 = smhaqsinha^, (30) 

and K' = i + K we obtain the final normalized operation 

£%p = e _x °{(c + eo + c-e 3 )p + (c+e 2 + c-ei)Z a pZ a 
+ (c+ei + c-e 2 )Z b pZ b + (c + e 3 + c^e )KpK 
+ (s+ei + s-e 2 )K'pK'i + ( s+e2 + s _ ei )K'^pK' 
+ (s+e 3 + s^e )JpJ J< + (s+e + S-e 3 )J^pJ}, (31) 



on the two-qubit state. Setting j3 b — [i a = [3 removes 
the single qubit terms Z a and Z b and also the K' terms, 
yielding the operation 



Up 



= e- 2 ^ 



(c+p + c-KpK + S-JpJ f + s + J f pJ) (32) 



with X3 — rift 2 . Focusing on the low loss regime, we can 
assume a small r\ between each interaction. Truncating 
to second order in r\ the expansion (26) leading to the 
above operation we have 



ri/3 2 (Z a + iZ b )p(Z a - iZ b ) 
r](3 2 Z a pZ a + rj0 2 Z b pZ b 



P 

irjP 2 (Z b pZ a - Z a pZ b ) 



(33) 



Let us assume we are in a quantum error correction 
(QEC) setting where ancilla qubits are being used. Once 
the ancilla systems have undergone projective measure- 
ments, syndromes are extracted for each logical qubit, in- 
dicating whether or not it has been subject to a Z error. 
The last two terms however will lead to cross terms of 
syndrome states which cannot be observed in the mea- 
surement process and thus they are removed from the 
resulting density matrix [56J. This leaves us with single 
qubit errors on a and b with equal probability (they con- 
stitute the observable part of the map). So for small loss 
we are only observing single qubit errors throughout the 
gate (at each one of the three dissipation stages) which 
can be corrected for via QEC. 

Correlated errors represented by the operator K ap- 
pearing in higher order terms are quantified by the coef- 
ficient c_ in the normalized operation (32). In Fig. 7(a) 
we can see how this part scales with loss in comparison to 
uncorrelated errors. The general quantum operations we 
have obtained are of great importance for error correc- 
tion. They provide us with the error syndromes and their 
associated probabilities which in turn can be directly fed 
into fault tolerance calculations. This will ultimately al- 
low us to compare the qubus scheme with other proposed 
implementations of quantum information processing. 

To appreciate the effect of loss on the whole gate we 
look at the fidelity of the output state with regards to the 
ideal output, for a two qubit input state |+)|+). Such 
an equally weighted superposition input state provides 
a good general indication as to the gates performance 
in addition to the fact that the ideal output is a useful 
resource, locally equivalent to a Bell state or two-qubit 
cluster state. A plot of the fidelity against the relative 
intensity decrease of the probe through dissipation de- 
fined as Itot = 1 — (exp[— 3^]) 2 is shown in Fig. 7(b). This 
can be understood in that if the probe is initiated in a 
coherent state with amplitude a, then at the end of the 
gate it disentangles from the qubits and is left with an 
amplitude ae~ 31 . In computing the fidelity we use 



Pa = Pb = r 





(a) 




0.07 






0.05 


S + + S. 
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FIG. 7: (a) The scaling of the normalized correlated (c_) 
and uncorrelated errors (s+ + s_) against the loss Itot as de- 
fined in the text. Dissipation occurs three times, once be- 
tween each interaction, in equal amounts, (b) The fidelity 
(dashed) F{\(f>),p) = {<j>\p\<j>) of the two qubit output state, 
where \4>) = e i7vZlZ2/4 \+} 1 \+) 2 constitutes the ideal output 
state and the concurrence (solid), both for the single and the 
iterated sequences. We observe the same fidelity for both of 
them while the outputted concurrence is lower for the iterated 
sequence. 



2 V e~ l + e 



-31 



(34) 



so as to make the phase represented by the last term in 
(22) that of an ideal CZ gate. We find that even for 
an intensity decrease of up to 80% in the probe mode 
(corresponding to 8 dB loss in total), the output fidelity 
remains above 0.5, allowing for purification. In the mod- 
erate loss regime, both the fidelity and the entanglement 
remain high; for example taking l to t = 0.05 correspond- 
ing to a decrease of 5% in probe intensity (0.22 dB loss) 
results in an output with F ~ 0.97 and a concurrence of 
-0.95. 



Finally we now show how a simple repetition scheme 
can significantly simplify full gate operation. The first 
point to notice is similar in spirit to the observation made 
in [39( that the ideal (loss-free) operation is invariant un- 
der time reversal. That is if we reverse the order of the 
interactions, we obtain the same conditional geometric 
phase. However in the case of separate cavities, a time 
reversed iteration does not help fight decoherence for a 
coherent state probe. The reason for this is that the sin- 
gle qubit error linked to the transfer between each cavity 
scales in the same way as the geometrical phase (of order 
2/3 2 ). The observation we make here is that the gate is 
also invariant under a swapping of the displacement di- 
rections. That is the same geometrical phase is obtained, 
again in the loss free case, if now qubit a conditionally 
displaces the probe in the imaginary direction in (22) and 
qubit b displaces it in the real direction. Let us denote 
the two different sequences in a dissipative setting as 



S = V{-i'l3 b Z b ).C.V(-'f3 a Z a ).C.V{i[3 b Z b ).C.V{(3 a Z a ) 

S = V(-i3 b Z b ).£.V(-i(3 a Z a ).£.V(p b Z b ).£.V(if] a Z a ). 

(35) 

We have 



S = e i < z ^-^' z ^ (66^)^(0), 



(36) 



with k — $ a Pb + PaPb- Now we notice that the effects 
of the two sequences differ only in the central overlap 
terms £2 and £2 which contain the multiple correlated er- 
ror terms. It is straightforward to see that their combined 
effect yields 



£26 = (i$ a Z a > + P b Z V \iJ3 a Za + fi b Z b ) 

x a z a > + i/3 b z b > \/3 a z a + if3 b z b ) 

= e ~^(i~^z b ,) e ~20l(l-z a z a ,) 



(37) 



which are just single qubit phase flip channels. Thus the 
combination of the two sequences gives the operation 



SS 



_ e 2iK,(z a z b -z a ,z b ,) 

x e -2vW b +l3 b )(l-^z„,) e -2 V (Pl+l3l)(l-z a z a ,) 



Z>(0), 



(38) 



where the first term is the unitary conditional two-qubit 
phase and the next two are single qubit dephasing terms. 
The corresponding operation undergone by the qubits, 
omitting the conditional phase is 



SSp 



(1 -Pa)(l -p b )p+p a (l -p b )Z a pZ a 

(1 - Pa)PbZ b pZ b + p a p b Z a Z b pZ b Z a , (39) 



with 



-*0(Pl\b+fi\ b ) 



Pa\b 



(40) 



the probability that each qubit incured a Z error. Clearly 
these dephasing processes are independent, leading to a 
very simple operation, in contrast with the partial oper- 
ations (31) or even (32) obtained in the single interaction 
sequence. Due to the fact that the geometrical phase and 
the qubit dephasing both double for a single iteration of 
the gate, there is no advantage to further reducing the 
sizes of the displacements and increasing the number of 
iterations. 



The output state fidelity and entanglement of this se- 
quence are plotted in Fig. 7(b), setting (3 a = (3 b = 



vvi 



e 3l ). The fidelity is very similar while the 



concurrence is slightly reduced, compared to the single 
sequence. So for entanglement distribution in view of 
communication applications, which only follow the pu- 
rification of an entangled state up to an acceptable level, 
the iterated scheme is penalizing, as it requires twice the 
amount of time and reduces the entanglement. But in 
view of full blown quantum computation with quantum 
error correction, the iterated scheme may present a seri- 
ous advantage, suppressing correlated errors at the gate 
level. 



V. DISCUSSION 

In this work we set out to examine the vital issue of 
dissipation in a CV quantum bus being used to real- 
ize quantum information processing tasks over discrete 
subsystems. We first looked at dissipation in the probe 
during interactions. We considered a general interaction 
Hamiltonian between continuous and discrete subsystems 
which appeared in two different physical settings: all op- 
tical cross-phase and dispersive CQED settings. Solving 
the master equation we derived the closed form of the co- 
herence parameter (7) discussed in some previous work 
[40|. This parameter can be expressed as a function of 
a, 7/x and \t and has a non-zero limiting behavior for 
large interaction times. As expected, it has the effect of a 
phase flip channel (10) in the case that the subsystem is a 
qubit. More interestingly, loss on the probe mode induces 
a conditional phase operation on the qubit, which needs 
to be undone in the case that the interaction is part of a 
more elaborate gate. The resulting time dependent den- 
sity matrix we were able to rewrite as a two-qubit density 
matrix. From this we observed the evolution of the en- 
tropy and the entanglement between the continuous vari- 
able and the qubit during the interaction. Increasing the 
probe amplitude results in higher entanglement with a 
sharper peaking, meaning the interaction time has to be 
controlled more and more precisely. We also found that 
increasing the probe amplitude increases the coherence 
of the combined state at entanglement peaking time. 

This could be a possible approach to minimizing deco- 
herence effects in qubus gates, however it has to be bal- 
anced against loss in between interactions and the con- 
ditional phase discussed above which depends strongly 
on the amplitude. As we saw in sections III and IV the 
dephasing effect on the qubit is related to the overlap 
between the possible states of the probe in phase space, 
which is a function of the squared distance between them. 
Thus increasing a in the interactions will augment the 
effects of loss in the transmission. The solution to this 
trade-off will depend on the physical system at hand. 

Having solved the effects of loss on interactions, we 
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combined them with displacements to simulate condi- 
tional displacements. These couplings are fundamental 
to hybrid schemes making use of the geometrical phase 
acquired by the probe in phase space. We managed to 
compensate for the combined dephasing effects during 
and in between interactions, obtaining the desired cou- 
pling up to a single qubit dephasing. This result is cru- 
cial to the analysis of the full CZ gate undertaken in 
section IV, in which we compounded the conditional dis- 
placement operations and added loss in between them 
too. The first and last loss terms constitute single qubit 
errors whereas the middle term contains correlated er- 
rors. We calculated the resulting quantum operation on 
the pair of qubits and found that in the low loss regime 
it only contained single qubit dephasing errors, given a 
QEC setting with syndrome measurements. To gain more 
insight into the gate's performance, we looked at the out- 
put state fidelity, for equally weighted input superposi- 
tion states. We found that the fidelity remained high 
and that even for losses of up to 8 dB in probe signal, 
the output states could be purified. A simple iteration 
procedure then allowed us to remove extra correlated er- 
rors altogether, presenting a serious advantage in view 
of a full-blown quantum computing scheme. This full- 
blown scheme will also have to factor other sources of 
errors such direct qubit decoherence processes. 

The method we use here is very broad and could be 
applied to general bus mediated quantum information 
processing schemes. Loss of coherence in the bus will au- 
tomatically result in dephasing on the qubits coupled to it 
at that time. If several qubits are simultaneously coupled 
then correlated errors will arise. These will have to be 
minimized in order to efficiently correct for errors. How- 
ever the peculiarity of the qubus scheme is that we are 
using non-orthogonal states of the bus to encode informa- 
tion held by the subsystems. Surprisingly, in the absence 
of loss, perfect multiqubit gates can be implemented this 
way. In a dissipative setting, this non-orthogonality leads 
to overlap calculations which are very likely to arise in 
general non-orthogonal bus schemes. 

Following the results presented in this paper, there are 
two main directions for future work. The first is to find a 
way of dealing with the single qubit phases accumulated 
throughout qubus gates. Simply undoing them would re- 
quire considerable precision and this may not be a realis- 
tic option. The second is to investigate possible schemes 
to reduce dephasing effects, possibly by engineering the 
probe itself. 



ACKNOWLEDGEMENTS 

We would like to thank S. Devitt, T.D. Ladd, P. van 
Loock, N. Liitkenhaus and T. Karasawa for useful in- 
sights and discussions. SGRL acknowledges the support 
of a Monbukagakusho scholarship. This work was sup- 



ported in part by MEXT in Japan and the EU project 
QAP. 

APPENDIX 

The exponent f nm in the coherence parameter ( nm = 
e /n m can kg wr itten as Re[/ nm ] + ilm[/ nm ]. Here we 
give expressions for the real and imaginary parts of the 
exponent with A = Z. Beginning with the real part we 
obtain 



,Re[/» 



exp[ 



;(X 2 (1 



,-27* 



2(7 2 + X 2 )' 
- 27 2 e- 27 *sin 2 x£ - X7e" 27 *sin2xi) 
x (l-Z n Z m )], (41) 

representing the decay of the off-diagonal components of 
the density matrix. In the limit of large interaction times 
this term leads to a fixed dephasing effect 



Re[/„, 



exp 



o ■ 



" \ J- %n ^m ) 



2(l + (7/x) 2 , 
Then the imaginary part of the exponent is 



(42) 



.M/] 



«7|a| 



eXP U( 2 i 2\ (*( 1 

2(7 +X) 
7sin2xi)(^ n - z m )], 



-27* 



cos2xt) 



(43) 



corresponding to a known single qubit phase term. 
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